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Abstract — We address the problem of a multi-user relay net- 
work, where multiple single-antenna node pairs want to exchange 
information by using a multiple antenna relay node. Due to the 
half-duplex constraint of the relay, the exchange of information 
takes place in two steps. In the first step, the nodes transmit 
their data to the relay, while in the second step, the relay is 
broadcasting the data by using linear and non-linear precoding 
strategies. We focus on the second step in this paper. We first 
consider the problem of maximizing the overall rate achievable 
using linear and dirty-paper type precoding strategies at the 
relay. Then, we consider minimizing the total power at the relay 
subject to individual SINR constraints using the same strategies 
at the relay. We show that the downlink-uplink duality does 
not hold for the setup considered here, which is a somewhat 
surprising result. We also show that the beamforming strategy 
which is optimal in the single-pair case performs very well in 
the multi-pair case for practically relevant SNR. The results are 
illustrated by numerical simulations. 

I. Introduction 

Cooperative communication plays a major role in future 
wireless ad-hoc -networks as well as in cellular systems such 
as LTE-Advanced. Cooperation can take place between base 
stations or mobile stations directly as well as via a relay 
station. Relay stations can hereby vary in the capability and 
complexity. They can be equivalent to base stations, i.e. with 
connection to the backbone, or as simple as a mobile station 
itself and variations between those extremes. The relays are 
then deployed within a cell or network in order to extend 
coverage or increase the battery life of mobile nodes. There 
is huge activity in the research community analyzing different 
aspects of relay networks. However, most of the work in multi- 
user relay networks is focused on the uni-directional case as 
e.g. in (JJ, (21 and is often limited to the case where multiple 
source nodes transmit data to their receiving counterpart(s) by 
exploiting the relay. 

Alternatively, the relays can also be used in order to 
exchange information between two nodes. The exchange of 
information between two nodes, often referred to as bidirec- 
tional communication, has been analyzed already by Shannon 
in |5). Some achievable rate regions for the bidirectional 
relay channel using different strategies at the relay, such as 
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decode-and-forward, compress-and-forward, and amplify-and- 
forward, have been analyzed in J4), Q. The performance of 
the bidirectional relay channel using superposition or network 
coding type strategies was analyzed in |6]]-|[8|]. In J9], some 
work has been done on near optimal relaying strategies as well 
as approximating the capacity region of the noisy (Gaussian) 
bidirectional relay channel. There, the relay uses a equation- 
forwarding scheme, in which the relay re-orders the received 
superposed signals, quantize and forwards them. It was shown 
that this scheme achieves rates which are within 3 bits of the 
cut-set upper bound. The scaling of the capacity region with 
multiple antennas and multiple relays for two-way relaying is 
considered in ITOl . 

For practical reasons, the relays involved in such cooperative 
communication systems often have an half-duplex constraint, 
i.e. they are not able to transmit and receive at the same time 
or at the same frequency simultaneously. Due to this half- 
duplex constraint of the relay, the exchange of information 
takes place in two steps. In the first multiple-access step, the 
nodes transmit their data to the relay, while in the second step, 
the relay is broadcasting data to the nodes. The capacity region 
of the broadcast step of this bidirectional relay channel was 
recently characterized in JT). 

Naturally, in a network there are multiple transceiver pairs 
which would like to exchange information by using the same 
relay. Obviously, this causes interference and thus in order 
to to have a reliable communication, interference manage- 
ment strategies are required IfTTl . IfTZl . For example, in IfTJl , 
ffl4l interference management is performed by allocating 
spread sequences and power efficiently to the communicating 
transceiver pairs. In |[T5l the capacity region of an determinis- 
tic lfT6l multi-pair bidirectional relay network is characterized 
and it was shown that the cut-set upper bound is tight. The 
results were then extended to the Gaussian case in |JT7| and 
an approximate capacity characterization is provided. 

The research on multi-pair bidirectional relaying is so far 
focused on single-antenna nodes. In this work, the relay is 
equipped with multiple antennas, while the nodes are equipped 
with a single antenna. With the generalization to a MIMO 
relay we will face new and interesting challenges, as we will 
see later on. Furthermore, the focus in this paper is on the 
broadcast phase of the bidirectional relay network. In more 
details, we are investigating the performance of two trans- 
mit strategies namely linear and nonlinear (dirty-paper |[T8l ) 
precoding performed by the relay in the broadcast phase. Our 
contribution is thus the performance analysis of these strategies 
in terms of achievable rate region using different beamformers: 



(a) First hop 



(b) Second hop. 



Fig. 1. Bidirectional relaying with multiple pairs. Pairs are indicated by 
equal numbers. 




(2.2) 



Fig. 2. Messages in the two-pair case 



• single-pair beamforming by treating the interference of 
the other (or remaining) pairs as noise 

• beamformers obtained after solving a related power min- 
imization problem with SINR constraints using a relax- 
ation method 

• an exhaustive search over a huge set of randomly gener- 
ated beamformers. 

Note that the downlink in the second phase is fundamentally 
different then the broadcast or multicast channel, since the 
receiving nodes have side information in form of their own 
messages conveyed in the first phase. 

The remainder of the paper is organized as follows. In sec- 
tionHU the system model is introduced, followed by section IITI1 
and [TV] in which the transmit strategies and the optimization 
problems are described. In addition, the lack of duality is 
discussed in this section. The results are then illustrated in 
section[V] followed by some concluding remarks in section [VTI 

II. System model 

Suppose that there is a group of single-antenna transceiver 
pairs communicating to each other in a bidirectional way by 
exploiting a relay equipped with N antennas as shown in 
Fig. 03 We assume that there are in total N pairs of nodes. 
Each node is denoted by the tuple (i, k), with 1 < i < N 
and k = {1, 2}, where the first index i identifies the pair and 
the second index k identifies which node of this node-pair is 
meant. Communication takes place in two hops. In the first 
hop, the multiple access step, the nodes transmit their data 
to the relay station. The relay station processes the received 
data and in the second hop, the broadcasting step, the data is 
forwarded to the nodes. 

In this paper, we analyze the broadcasting step by using 
linear and nonlinear precoding strategies at the relay. For 
simplicity, we assume that the first hop was successful, i.e. 
all messages were received with an asymptotically small error 
probability. This assumption is valid given the rates are within 
the MAC capacity region in the first hop and the first hop does 
not pose a bottleneck for the system. The resulting commu- 
nication scenario is illustrated in Fig. [2] for the two-pair case. 
In this case, the relay has in total four independent messages, 
given by W^i), W^(i,2)> W(2,i)> an d W( 2 ,2)- Obviously, each 
node is aware of its own message which was transmitted in the 
first hop. Furthermore, each node is interested in the message 
from the other node belonging to the same pair, illustrated by 



From the information available at the receivers, it becomes 
clear that this problem distinguishes itself from a multicast 
setup, in which a common information is broadcasted to a 
group of receivers without a priori information at the receivers. 
For the same reason, it is also different from the well studied 
broadcast channel. We would also like to emphasize that there 
are in total TV pairs of nodes, i.e. we have 2N nodes in 
total, while the relay is only equipped with N antennas. As 
a consequence, complete interference avoidance using zero- 
forcing transmission is not possible. 

In the following, we are dealing basically with two op- 
timization problems. In the first optimization problem, we 
are optimizing a weighted sum-rate subject to a transmit 
power constraint at the relay. In the second optimization 
problem, we minimize the total transmitted power from the 
relay under a minimum SINR requirement for each node. For 
both optimization problems we consider a linear as well as a 
non-linear precoding strategy. 

III. Maximizing sum rate 
We start with the linear strategy. 



A. Linear techniques 

Using linear precoding at the relay, the received signal at 
node (i, k) is given by 

V{i,k) = V^ h (a-) u ' s * + E \ ','>. + n (i,k), (1) 
for« = l,...,iV, k = l,2, (2) 

where pi is the power allocation to node-pair i, Ui is the unit- 
norm beamformer for node-pair i, h(,-M describes the channel 
to the fc-th node, k = 1,2, of node-pair i. The channels 
between the relay and the nodes are modeled as complex 
Gaussian random variables, i.e. CAf(Q, 1). The relay has an 
power constraint Y^,iPi < P, n>ujA is the additive white 
Gaussian noise with variance a 2 at the fc-th node of node- 
pair i. Thus, the transmit signal is of the form 
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where s, is the complex valued signal denoting the information 
signal for node-pair i. The received SI N fe -, at node (i, k) is 

2 

Pi h? 
SlNR (iifc) 



(4) 



It has been shown in (6), Ifl9l for the special case of one pair, 
that the capacity region can be characterized by its boundary 
which corresponds to the weighted sum rate given by 



R = M(i,i) lo S2C 



M(l,2 



C 
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(5) 

with weights /Un i) and /X(i,2)> where C(x) = log 2 (l + x). 

Thus, using Q and treating the multi-user interference as 
noise, the sum rate can be characterized for the general case 
by 



i? = ^C(SINR (l , 1) )+C(SINR (ij2) 



(6) 



i=l 



For the two pair case the sum rates of each pair are given by 



Ri 




R 2 = C 



In the following paragraph, we are deriving the rate expres- 
sions achievable with non-linear precoding, namely dirty paper 
coding. We then describe the choice of beamforming using 
non-linear precoding with the goal of maximizing the weighted 
sum rate subject to power constraints. The differences in the 
expressions for linear precoding are described along with the 
discussion on the algorithm. 

B. Dirty-Paper Techniques 

With dirty -paper coding |fl~8), we are able to pre-compensate 
the interference since it is already known at the transmitter. By 
assuming an arbitrary but fixed encoding order, the SIIMR^) f. 
of user k of pair i is given by 
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h (l(i),fc) U 7r(j) 



(7) 



Thus, we have the following rate vector R = 



[-R(E,7T( 



1)), • ■ • ! R(S,-K(i)) , 



; -R(e.7t(jv))]i where the individual 

(8) 



sum rates i?(s iJr (i)) for each pair are given by 

R(s,n{i)) = C (SINR( 7r (i) j)) + C (SINR (7r( j)^) 
The dirty rate paper region i?opc is then given by 

-Rdpc = ConvexHull I [J R(7T,pj) J , 



(9) 



which is the the convex hull of the union of all rate vec- 
tors R(7r,Pi) over all powers pi and over all permutations 

(7r(l),...,7T(JV)). 

Assuming that the interference of the other nodes in (O 
and (O can be regarded as Gaussian, in our approach the 
beamforming vectors are computed very efficiently for a given 
order of precoding in the following successive way. First of 
all, from expression (0, we observe that the node pair n(N) 
encoded last experiences no interference from the other nodes. 

Thus, using the approach in |fl9l the beamforming vector 
for the node pair tt(N) is obtained by 

M tSU(N).l) + (!-*) CXp(-J0)g( ff (iV),2) 

u ir(JV)W — Tu m 7\ 7 — 7\ TT { - lv > 

ll*g(7r(jv),i) + (1 ~ t)exp(-j(f))g i7T (N),2)\\ 



l (Tr(N),l)/\\h(Tr(N),l)\\' S(ir(iV),2) 



with g( OT (jv),i) 

\n(N),2)/\\h(„(N),2)\\, and (f) = arg ^g^. (JV)jl) g (w(JV)j2 
The weight t can be used to prioritize one of the nodes 
Intuitively, the beamforming vector in dTOb represents a linear 
combination (up to a phase adjustment) of maximum-ratio 
transmission beamformers in the direction of the node pairs. 

The node pair tt(N — 1) is considered next. The node pair 
ir(N— 1) observes interference from the beam intended for the 
node pair tt(N), which is exactly known. Assuming again that 
the interference from the node-pair tt(N) can be regarded as 
additional white Gaussian noise then the overall interference 
plus noise (referred to as effective noise in the following) is 
distributed as CJ\f(0, er? 



(7T(N-l),fc) 



with 



r (7r(iV-l),fc) — Ptt{N) 



n (7r(AT-l),fe) U 7r(A r ) 
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we compute the beamforming vector for node pair ir(K — 1) 
very efficiently as in ( TTOb . 

This process is continued until is determined. The 

above algorithm has to be repeated for all possible user 
orderings it. In the linear precoding case, each receiver is 
observing interference from all other beams in the system, 
thus the effective noise variance (fTTT i is changed accordingly. 
The beamformers are obtained by evaluating (TTOb accordingly. 
Note that there is no claim of optimality in terms of achievable 
rates of the beamforming and precoding approaches discussed 
above. Due to the non-convex structure of the optimization 
problem, an optimal beamforming strategy for the weighted 
rate maximization is difficult to obtain. However, note that 
beamforming has always the advantage of simple processing 
at the transmitter (scalar instead of vector coding) and receiver 
(single stream decoding) and is thus often considered in 
wireless standards such as LTE and WiMAX, which justifies 
the analysis of such schemes. As an alternative to the ap- 
proach discussed above, in the following section the equivalent 
problem of power minimization subject to quality-service- 
constraints (QoS) in terms of required SINR is considered, 
first for the linear precoding case, followed by non-linear 
precoding. 



IV. Minimizing power 

A. Linear techniques 

In the following, we will show that the uplink-downlink 
duality lEOl . ED does not hold in the setup considered here. 
We start with the downlink. The optimization problem by 
using linear precoding is given by 

N 



(12) 



s.t. 3ui,...,ujv,||ui 



1,1 < I < N, 
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7(i,fe), 



k = l,2,i = 1,...,N 



where 7(j,fc) is the SINR requirements at node (i,k). Let us 
introduce the following matrices with 



V (fc) _ \ h H . 



V (fe) = 0, 



1 < i,j < N,k = 1,2 (13) 
1 < i < AT. (14) 



Note that the V( fe ' are functions of the beamforming vectors, 
i.e. V( fe ) = V< fc ) (U), with U = [ui, . . . , ujv]. Then <0 can 
be rewritten in the following form 



N 

min. ^2 Pi 



(15) 



s.t. 3U, \\u,\\ 2 = 1,1<1<N 

P > (p^y 1 (v^V^p + a^h) k = l,2 
D (fc) = diag 



with r' fe ' = diag (7(1, fc), • ■ ■ , 7(iv,fc)) representing the matrix 
of SINR requirements, where 1 is the all one vector. The power 
minimization problem in the dual uplink is given by 
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■E' 

i=l 



(16) 



t.q> D« rW V") q + ^l fc = l,2 



Note that at the optimum the constraints are achieved with 
equality (otherwise, the power corresponding to the constraint 
f° r 7(i,fc) can t> e reduced such that equality is obtained, 
which would reduce the power spent and thus would be a 
contradiction.) In the following, by using a counter-example 
we show that uplink-downlink-duality does not hold. Consider 
the case N = 2, with = I, k = 1,2, = I, and 

r^ 2 ) ^ I. Then, define the following sets according to the 
constraints in (TT3T > and ( TTST l for this case 



-Mdl = {p G E.+ : p > r (fc) V (fc) p + cr 2 r« [ I] , k = i, 2} 



(17) 



Mul = {q G R+ : q > T (fc) (V (fc) ) H q + a 2 T (k) [ \ ] , k = 1, 2} 




1 2 3 4 Pi 1 2 3 4 9i 

(a) Downlink (b) Dual Uplink 

Fig. 3. Illustration of counter-example to uplink-downlink duality 



for downlink and uplink, respectively. Note that .Mdl = 
Hfe-M^L an( l similarly for the uplink. The sets are illustrated 
in Fig. |3 In Fig. 0a) (analogously in Fig. Ob)) there are 
three relevant intersections, marked with a number (1), (2), 
and (3) respectively. The numbers (1) and (2) indicate the 
minimum of the subset and M^, respectively. The point 
(3) indicates the minimum of the Mdl- For simplicity, it is 
assumed that V"W is symmetric (cf. Fig. [3}. It follows that 
A^e>l ' s s y mmeu "i c an( l so is -^{jl- Note that the minimum 
of .Mql an£ l -^ul ^ on tne same dotted slope —1 line 
passing through (1) due to duality (since £ 4 Pi = £j ?«)■ The 
same holds for M^l and -M{j L . However, as can be seen in 
the figure it holds not true for .Mdl a nd .Mul (compare the 
dotted lines passing through the points (3)), unless = I 
as well, which is in general not the case. 

Thus, uplink-downlink duality can not be exploited in order 
to get a solution for the downlink beamformers as in ll20l . OTl 
for the broadcast channel (unless the SINR requirements r( fe ) 
are equal for all nodes, which is a case of limited interest). 

However, the optimization problem in (fT2l can be solved 
approximately by using a semidefinite relaxation. Let us define 
Qi = piUiuf . In the original problem, the Q^, Vi, are rank- 
constraint. This non-convex constraint can be replaced by a 
convex constraint >z 0, where the notation A >z B means 
A — B is a positive semidefinite matrix. While the original 
constraint was of rank 1, the Qi can be of any rank with 
this relaxation. Thus, the relaxed problem is a semidefinite 
programming problem (SDP) given by 



TV 



minimize trace (Qj) 



(18) 



i=l 



subject to 



_ 7(i,fc) (19) 



> 



t 0. (20) 
The solution to the relaxed SDP (fL8b gives a lower bound of 



the objective function and a relaxed solution Q*. However, if 
that solution Q* is of rank 1, then the solution to the relaxed 
optimization problem is identical to the solution of the original 
problem. In general, however, the solution will not have rank 1. 
As an heuristic but efficient approach, here the eigenvectors 
corresponding to the largest eigenvalue of Q* are used as a 
rank 1 approximation. 

B. Dirty-paper techniques 

The procedure here is similar to the linear beamforming 
case, except that the order of encoding is an additional parame- 
ter which has to be taken into account in the optimization. The 
SDP relaxation in the case of non-linear precoding is given by 

N 

minimize trace (Qi) (21) 
i=i 

,. tt h Mi),fc)QT(i) h (>(i),fc) 

sub J ect to v hH n h T^2 ^ 7«i),fc) 

(22) 

Qnii) h 0. (23) 

Once the optimal beamformers are obtained from i2li . a 
rank 1 approximation is obtained using the same procedure as 
in the linear case. By defining the SINR requirements Y^jVfc) 
in terms of rates as follows 

70r(i),k) = 2^W fl =-"W - 1 (24) 

with J^i Mi = 1> we can a PPly an iterative (bisection) El 
power minimization problem to get an achievable rate region. 
The parameter Rs, v U) is changed at each iteration, i.e. the 
value of -Rs.Tr(i) is incremented in the next iteration if (fJTJ 
is feasible otherwise decreased. The step size is obtained by 
bisection between the last feasible value of Rs^U), which is 
initially equal to zero, and the last infeasible value of R^^u), 
which can be chosen as twice the capacity in the single-pair 
case given in (O for initialization. The iterations continue until 
a predefined accuracy is achieved. 

V. Illustration 

In this section, the results are illustrated by means of 
numerical simulations. Due to the similarity of the results, 
we focus in this section on the non-linear precoding strategy, 
i.e. dirty-paper coding. In Fig. @] the rate region of a two 
pair network is shown for a SNR = 3 dB. Note that on 
the coordinate we have the sum rate of pair i = 1, while 
on the ordinate we have the sum rate of pair i = 2. Thus, 
the figure is a projection of a four-dimensional rate region 
to a two-dimensional one. For both curves in the plot, dirty- 
paper coding is used. As beamforming strategy the single- 
pair beamforming approach (cf. ([Tol l) is used and compared 
to an exhaustive search over a huge set of randomly generated 
beams. The channel vectors used in the plots have the follow- 
ing norms: ||hu||| = 3.28, ||hi 2 ||| = 2.9,||h 2 i||| = 1.77, 
and HI122II2 = 2.2. From the figure, we observe that the single- 
pair approach performs very close to the exhaustive search. 




Fig. 4. Rate Regions achieved with Dirty-Paper Coding; single-user beam- 
forming approach (cf. (To)) vs. randomly generated beams at SNR = 3 dB 




Fig. 5. Rate Regions ach ieved with Dirty-Paper Coding; single-user beam- 
forming approach (cf. i\0\ ) vs. randomly generated beams at SNR = 10 dB 



The same strategies are compared in Fig. [5] now for 
10 dB. Here, we observe a small gap between the single- 
pair approach and the exhaustive search. Still, the single-user 
strategy performs reasonably well. 

The next Fig. [6] shows the performance of the two strategies 
for 30 dB. In addition to this, the region obtained using the 
relaxation method given in $2l[ is plotted. From the figures, 
we observe that the gap between the exhaustive search and the 
single-pair approach has increased further. The region obtained 
using the relaxation method is comparable to the one using the 
single-pair approach. From the figures, we conclude that our 
approach is close to the exhaustive search for low and average 
SNR, but is suboptimal for high SNR. 

VI. Conclusion 

In this paper, we have addressed the problem of a multi- 
user relay network, where multiple single-antenna node pairs 
want to exchange information by using a multiple antenna 
relay node. In the first step, the nodes transmit their data to 
the relay, while in the second step, the relay is broadcasting 
the data by using linear and non-linear precoding strategies. 
We focused on the broadcasting step and first considered the 
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Fig. 6. Rate Regions achieved with Dirty-Paper Coding; single-user 
beamforming approach (cf. fllOt ), randomly generated beams, and the re- 
gion obtained (dashed-dotted line) using the relaxation method (cf. )2U ) at 
SNR = 30 dB 



problem of maximizing the overall rate achievable using linear 
and dirty-paper type precoding strategies at the relay. Then, 
we considered the problem of minimizing the total power at 
the relay subject to individual SINR constraints. We showed 
that the downlink-uplink duality does not hold for the setup 
considered here. We also showed that using the beamforming 
strategy which is optimal in the single-pair case performs very 
well for practically relevant values of SNR. 
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